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ABSTRACT 

We analyse absorption line-strength indices for ^3000 red-sequence galaxies in 94 
nearby clusters, to investigate systematic variations of their stellar content with lo- 
cation in the host cluster. The data are drawn from the NOAO Fundamental Plane 
Survey. Our adopted method is a generalization of that introduced by Nelan et al. to 
determine the global age-mass and metallicity-mass relations from the same survey. 
We find strong evidence for a change in galaxy properties, at fixed mass, over a range 
from the cluster centre to the virial radius, i?200- For example, red-sequence galaxies 
further out in the clusters have weaker Mgb5177 (at ~8c significance) and stronger 
H/3 and H7 absorption (^3<r, ~4cr) than galaxies of the same velocity dispersion in the 
cluster cores. The Fe5270 and Fe5335 indices show only very weak trends with radius. 
Using a total of twelve indices, the pattern of cluster-centric gradients is considered 
in light of their different dependences on stellar age and chemical composition. The 
measured gradients for all twelve indices can be reproduced by a model in which red- 
sequence galaxies at ~li?200 have on average younger ages (by 15 ± 4%) and lower 
cv-element abundance ratios (by 10 ± 2%), than galaxies of the same velocity disper- 
sion but located near the cluster centres. For the total metallicity, Z/H, no significant 
gradient is found (2 ±3% larger at -R200 than in the cores). There are hints that the age 
trend may be stronger for galaxies of lower mass and/or for galaxies with more disky 
morphology. We show, however, that the trends cannot be driven primarily by changes 
in the morphological mix as a function of radius. The cluster-centric age and [a/Fe] 
gradients are in the sense expected if galaxies in the cluster core were accreted at an 
earlier epoch than those at larger radii, and if this earlier accretion contributed to an 
earlier cessation of star formation. The size of the observed age trend is comparable 
to predictions from semi-analytic models of hierarchical galaxy formation. 

Key words: galaxies: evolution — galaxies: elliptical and lenticular, cD — galaxies: 
clusters: general 



1 INTRODUCTION 

Evidence that galaxy properties depend on their large-scale 
environment dates from at least seven decades ago (Hubble 
& Humason 1931), with the morphology-density relation- 
ship being the classic manifestation (e.g. Melnick & Sar- 
gent 1977; Dressier 1980). The spectral properties of the 
galaxy population similarly exhibit environmental correla- 



tions, with more emission-dominated objects in the field 
than in clusters. Measurements of Ha equivalent widths in 
large spectroscopic surveys reveal a suppression of the aver- 
age star-formation rate for galaxies in clusters, with this ef- 
fect extending well beyond the virialized cluster core (Lewis 
et al. 2002; Gomez et al. 2003). A distinct but related issue is 
whether the correlation of star-formation rate with density 
is simply a consequence of the morphology-density relation, 
or whether there is a suppression for cluster galaxies relative 
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to field objects of the same morphological type (e.g. Balogh 
et al. 1998). 

A striking feature of the galaxy population within rich 
clusters is the sharp envelope in colour-magnitude space de- 
fined by the reddest galaxies, and interpreted as the locus 
of non-star-forming galaxies. Objects on this 'red sequence' 
are predominantly of early-type morphology, i.e. ellipticals 
and SOs; their structural and dynamical parameters follow 
tight scaling relations such as the Fundamental Plane (e.g. 
Djorgovski & Davis 1987), while their colours and spectro- 
scopic line strengths reveal that their average stellar popu- 
lations are correlated with the luminosity or mass (e.g. Vis- 
vanathan & Sandage 1977). After star formation ceases, a 
galaxy evolves onto the red sequence within a fairly short 
time interval of ~lGyr. This rapid reddening helps to ex- 
plain the strong bimodality of colours in the overall galaxy 
population (e.g. Hogg et al. 2002). 

Beyond one Gyr after the end of star formation, the 
broadband optical colours become less sensitive to the stel- 
lar age; moreover, the effects of increased metallicity closely 
mimic the effects of greater age. This classic 'age-metallicity 
degeneracy' allows the red sequence to conceal a wide range 
of ages, especially if there is a 'conspiracy' such that galaxies 
with younger stars also have higher metallicities, at given 
mass. To resolve age differences within the red sequence, 
indicators which are more sensitive and more robust than 
broadband colours are needed. A widely-applied technique 
exploits the age dependence of the Balmer absorption lines, 
in combination with metallicity-driven indices, to break the 
age-metallicity degeneracy (e.g. Worthey 1994). A number 
of recent works with this method have concluded that age 
must vary along the red sequence, with less massive galaxies 
having younger ages than more massive systems, (e.g. Cald- 
well et al. 2003; Nelan et al. 2005; Bernardi et al. 2006). 
These studies also support an increase of overall metallic- 
ity with mass, as well as an increase in the abundance of a 
elements relative to iron. 

Beyond the existence of the overall age-mass trend, the 
distribution of red-sequence galaxy ages should reflect how 
star formation ended in galaxies in different environments. 
In N-body simulations of hierarchical clustering, the average 
accretion time of sub-haloes in cluster-sized parent haloes 
appears to be correlated with their final radius from the 
cluster centre (Gao et al. 2004, but see also Moore, Diemand 
& Stadel 2004). In the simulations analysed by Gao et al., 
galaxies within 20% of the virial radius at z — entered 
the cluster at an average of z ~ 1 (lookback time ~8Gyr), 
while those at the virial radius were accreted only at z ~ 0.4 
(~4Gyr ago). Any mechanism which suppresses star for- 
mation when a galaxy enters a cluster, combined with the 
accretion-time gradient, will yield a cluster-centric gradient 
in the stellar populations. As discussed by Balogh, Navarro 
& Morris (2000), an attractive candidate for this mechanism 
is 'suffocation': the removal of the hot gas reservoir from the 
in-falling galaxy, causing a gradual decline in star-formation 
rate on time-scales of 1-3 Gyr (see also Larson, Tinsley & 
Caldwell 1980). If such a process acted on at least some of 
the galaxies now on the red sequence, a consequence will be 
for objects now at larger radius to have experienced more 
extended star-formation histories on average. At any given 
mass, the outer galaxies will have luminosity- weighted ages 
slightly younger, and a-element abundances (driven by the 



ratio of prompt SNe II versus delayed SNe la) slightly lower 
than galaxies near the core of the cluster. This behaviour 
is now captured qualitatively by semi-analytic models of 
galaxy formation, which combine numerical prescriptions for 
star formation, and other baryonic physics, with N-body 
dark-matter merger trees. For instance, in the models of 
De Lucia et al. (2006) , the mean luminosity- weighted age of 
bulge-dominated galaxies falls from ~ 12 Gyr at cluster cen- 
tres to ~ 10 Gyr at the virial radius. Although some of this 
effect is due to mass segregation (since the largest galaxies 
are slightly older, and preferentially located near the centre), 
this is not the dominant cause, and a true residual trend of 
age with radius, at fixed mass, appears to be present. 

While there is substantial observational evidence for 
truncation of star formation in clusters for the galaxy 
population at large, there is less conclusive evidence for 
environment-driven effects in the red galaxies themselves. 
In this case, since the mass-sequence of galaxies dominates 
most of the observed correlations, any environmental ef- 
fect would be revealed by systematic shifts in the residuals 
from the overall scaling relations, in particular the colour- 
magnitude relations (e.g. Terlevich, Caldwell & Bower 2001; 
Hogg et al. 2004) or the spectroscopic line index versus ve- 
locity dispersion (index— a) relations. Comparing samples 
drawn from field and cluster environments, Bernardi et al. 
(1998) found marginally smaller Mgb5177, at given velocity 
dispersion, for their field sample. Kuntschner et al. (2002) 
used a range of indices to infer that very isolated early- 
type galaxies have younger stellar ages than their counter- 
parts in clusters. An alternative to the cluster- versus-field 
approach is to compare the properties of cluster galaxies 
at different radii from the cluster centre. An early study of 
this kind was made by Guzman et al. (1992), who found 
a significant offset in the magnesium index Mg2, between 
the core of Coma and galaxies at ~ 1.5i?200. The results 
were interpreted as evidence for a gradient in stellar ages, 
with more recent star- formation in the halo sample. Also in 
Coma, Carter et al. (2002) have claimed to find an offset 
in the same sense for Mgb5177, but only very weak gradi- 
ents in the age-sensitive H/3 and the iron-abundance indi- 
cator (Fe) = i(Fe5270 + Fe5335). Considering the relative 
strength of the effect in the three indices, Carter et al. favour 
an explanation in terms of a cluster-centric metallicity gra- 
dient, and invoke pressure confinement of supernova ejecta 
by intra-cluster gas as the underlying mechanism. 

In this paper we determine the cluster-centric gradients 
in the index— a residuals, using data from the National Op- 
tical Astronomy Observatory (NOAO) Fundamental Plane 
Survey (NFPS). We use the gradients for a set of twelve 
spectral indices, with different responses to the stellar pop- 
ulation parameters, to make a clear separation of age and 
metallicity effects. The method applied is a generalisation of 
that used by Nelan et al. (2005) to investigate the trends of 
stellar populations with increasing mass, based on the same 
sample. The layout of the present paper is as follows: In 
Section [5] we first describe the sample and data employed, 
while in Section |3] we determine the radial trends followed 
by twelve individual line-strength indices. These results are 
independent of any adopted stellar population models. In 
Section^] we use models of Thomas et al. (2003, 2004) with 
the aim of describing the radial behaviour of all twelve of 
the indices, in terms of underlying trends in age, metallicity 
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and a-element enhancement. We consider in some detail the 
degree to which a morphology-radius trend contributes to 
these results, and present extensive tests for systematic bi- 
ases. Section [I>] compares the results obtained against previ- 
ous work within clusters, and with cluster- versus- field stud- 
ies. We also consider the physical origin of cluster-centric 
age gradients in hierarchical structure formation models, 
and compare our results with predictions from semi-analytic 
modelling. Finally, our principal conclusions are reiterated 
in Section [(J 



2 DATA 

The data to be used in this analysis are fully described in 
Smith et al. (2004), which also presents an overview of the 
NFPS programme, and in Nelan et al. (2005), where the line- 
strength measurements are reported in detail. We refer the 
reader to these papers for a full discussion; here we provide 
only a summary. 

The measurements are based on multi-fibre spectro- 
scopic data obtained at the Cerro-Tololo Inter-American 
Observatory (CTIO) 4-m telescope and the Wisconsin- 
Indiana- Yale-NOAO (WIYN) 3.5-m telescope. Absorption 
line-strengths, as defined by Trager et al. (1998) and 
Worthey & Ottoviani (1997), were measured after smooth- 
ing the spectra to match the resolution of the Lick stellar li- 
brary, in order to match predictions based on the Lick fitting 
functions. Unlike some studies, we made no attempt to bring 
the data onto the Lick system by comparing standard stars. 
Rather, the indices are measured on flux-calibrated spectra. 
As will be seen, our method employs only relative changes 
in the line-strengths, so absolute zero-point shifts from the 
Lick system would not affect our results. We do not in this 
paper apply aperture corrections to the line-strengths based 
on galaxy distance, preferring to absorb such effects into a 
free parameter in our fits. For the velocity dispersions, we 
apply an aperture correction of 0.04 dex per decade in angu- 
lar diameter distance (e.g. J0rgensen, Franx & Kjaergaard 
1995). 

The set of twelve indices analysed here is the same as 
used by Nelan et al. to determine the mass scaling rela- 
tions. As in the previous paper, we will interpret the index 
data in the context of models by D. Thomas and collabo- 
rators (Thomas, Maraston & Bender 2003; Thomas, Maras- 
ton & Korn 2004; jointly hereafter TMBK). These models 
are widely used since they incorporate the effects of abun- 
dance ratio variations, through the parameter [a/Fe]. The 
super-solar [a/Fe] measured in giant elliptical galaxies (e.g. 
Worthey, Faber & Gonzalez 1992) are interpreted as evi- 
dence for star-formation time-scales which are short com- 
pared to the ~1 Gyr lifetimes of Fe-producing SNIa. More 
recently, abundance ratio effects have been incorporated 
into alternative models (e.g. Lee & Worthey 2005; Schi- 
avon 2005); for ease of comparison, however, we use only 
the TMBK models in this paper. There is an element of 
subjectivity in the choice of indices to incorporate in the 
fits, but the index set needs to be sensitive to the three 
parameters covered by the models, i.e. age, metallicity and 
Q-abundance. Our chosen set includes three a-element in- 



dicators (Mgb5177, CN1, F64668 1 ), and three age-sensitive 
Balmer indices (HdF, HgF, Hbeta) 2 . The remaining six in- 
dices measure primarily the iron-peak metallicity (Fe4531, 
Fe5270, Fe5335, Fe5406, Fe5709, Fe5782). The balance be- 
tween the three groups is intended to compensate for the 
intrinsic sensitivity of the indices: the Fe indices only de- 
pend weakly on the underlying parameters, relative to the 
typical errors, so a greater number should be included to 
reach a given significance level. 

Because the NFPS galaxy sample covers a broad red- 
shift range, 0.01 < z < 0.08, the strong night-sky lines at 
5577 A and 5895 A preclude accurate line-strength measure- 
ments for indices redder than Mgb5177, over certain red- 
shift intervals. For instance Fe5270 cannot be measured for 
cz between 14500 km s -1 and 20000 km s~ . A feature of our 
analysis method, described below, is that we do not require 
an identical set of indices to be measured for all galaxies in 
the sample. Using multiple Fe line-strengths in particular al- 
lows the whole redshift range to be covered by multiple semi- 
redundant indices. In the case of the reddest indices, Fe5709 
and Fe5782, we further restrict the redshift range such that 
measurements at observed wavelengths longer than 6050 A 
do not enter the analysis. This region corresponds to the 
extreme red part of the spectrum in the CTIO data, and 
may be subject to flux-calibration problems. Note that for 
consistency between northern and southern data subsets, we 
reject all line-strengths beyond this cut-off, including those 
measured from WIYN spectra. 

Although the line-strength indices are designed to mea- 
sure stellar absorption features in integrated spectra, some 
of the indices are contaminated by the presence of nebular 
emission lines at the same wavelengths. The most trouble- 
some cases are for the age-sensitive Balmer indices, contam- 
inated by nebular emission at H/3, H7 and H<5. The effect 
of emission is to 'fill in' the stellar absorption, leading to 
erroneous age estimates which are older than would be de- 
termined in the absence of emission. Relative to the age- 
effect in the stellar spectrum, the nebular contamination is 
more severe for Ha and H/3 than for the higher-order lines. 
With data of very high quality this can in principle be made 
the basis of an emission-correction scheme (e.g. Caldwell, 
Rose & Concannon 2003). For NFPS, however, only a small 
subset of galaxies have Ha measurements (see Smith 2005). 
More traditionally, emission corrections have assumed a con- 
stant ratio of [O ill] to H /3 equivalent widths, and used the 
easily-measured [O ill] line at 5007A to estimate the contam- 
ination at H/3 (e.g. Gonzalez 1993). In Nelan et al. (2005) we 
showed, however, that the [O ill] and H /3 data for the NFPS 
sample show a wide (factor of ~20) variation in line ratios. 
Moreover the distribution of line ratios is bimodal (see also 
Yan et al. 2005), and strongly related to the galaxy mass 
or luminosity (as also discussed by Goudfrooij 1999). As 
a result, the frequently-adopted factor of 0.6x[Oiii A5007] 
(Trager et al. 2000) is a poor and misleading correction, 



1 Sometimes called C24668 because of its strong dependence on 
the C abundance, see Tripicco & Bell (1995). 

2 We will adopt the convention that HdF etc refer to the line- 
strength indices, defined by fairly broad continuum and central 
pass-bands, while H<5 etc refer to physically meaningful spectral 
lines, either in emission or absorption. 
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even in a statistical sense. In the present paper we apply 
a more conservative approach, selecting the galaxy sample 
to exclude objects with nebular emission according to the 
measurements of Nelan et al. Specifically, we adopt selection 
limits of EW(B./3) > -0.6A and EW [Om A5007]> -0.8 A 
(where the equivalent widths EW are negative for emis- 
sion, as tabulated by Nelan et al.), and galaxies not meeting 
both of these criteria are removed from the sample entirely. 
Thus all of our results will strictly refer to emission-free red- 
sequence galaxies. 

For an explicitly cluster-based survey such as NFPS, 
the projected cluster-centric distance, R c i, is a natural and 
convenient (although imperfect) parameter to describe the 
environment of a galaxy. We adopt the nominal cluster cen- 
tres as reported by Smith et al. (2004), which were based 
on either the position of one or more dominant galaxies, 
or a clear X-ray emission peak, on a case-by-case basis. To 
account for the range in richness of the NFPS clusters, we 
express R c \ in terms of the cluster virial radius, f?200, defined 
as 



-R200 = ( 



(Tel 



1000 kms 



h' 1 Mpc 



(1) 



where cr c i is the cluster velocity dispersion (e.g. Carlberg 
et al. 1997) and h is the Hubble constant in units of 
100 kms -1 Mpc -1 . We note that alternative definitions have 
been used for the virial radius by other authors (e.g. Quin- 
tero et al. 2005), but in the present paper we are concerned 
primarily with placing any gradients on a relative scale, and 
to first order the above should be sufficient. 

Finally, since this work will investigate radial trends in 
the properties of 'normal' red-sequence galaxies, we need to 
be careful in our treatment of cD galaxies. Since cDs are 
found preferentially near the cluster centre but probably 
have distinctive star-formation histories relative to other gi- 
ant elliptical galaxies (e.g. Crawford et al. 1999) , they could 
exert an undue influence on our fits. However, some cDs 
in our sample are not coincident with the nominal cluster 
centres as listed by Smith et al. (2004), so an inner radial 
cut is not sufficient to remove all of them from the sample. 
Instead, we use the J-band luminosity from the Extended 
Source Catalogue of the Two Micron All Sky Survey (Jar- 
ret et al. 2000), to exclude from the fit ~200 galaxies with 
Mj < —23.25 (for h — 1). For the isophotal magnitudes used 
for the selection, this limit is approximately one magnitude 
brighter than M} (Kochanek et al. 2001; Cole et al. 2001). 
Although the cut will remove some normal giant ellipticals 
from the sample, as well as the cDs, this should not bias 
the fits; the sample size is reduced by only <10%, causing 
minimal loss in the overall signal-to-noise ratio. 



3 MODELLING THE INDEX er RELATIONS 

In this section, we fit linear multivariate models for the index 
measurements, including a term describing cluster-centric 
variations at given mass. The results of this section are based 
only on observed variations in the line-strengths and hence 
do not depend on any particular stellar population model 
used to interpret the index data. 



3.1 Motivation 

Before presenting the fits in detail, we first illustrate the 
presence of cluster-centric gradients using three example in- 
dices, Mgb5177, HgF and Fe5270. Figure [I]shows the index- 
cr data for these indices, with the global relations as derived 
in the following section. As shown in Nelan et al., the rel- 
ative slopes of the index-cr relations can be understood in 
terms of systematic variations in stellar age and chemical 
composition along the mass sequence. 

As initial evidence for cluster-centric variations, Fig- 
ure|5|shows quantiles of the index— a relations in two ranges 
of Rci, corresponding to the cluster core (R c \/R200 < 0.2), 
and to the outermost galaxies observed (0.5 < Rd/Rioo < 
1.0). The global relations are reproduced in each panel for 
comparison. The shift in the Mgb5177— a relation is readily 
seen, and is well described by a simple offset at all a, with 
the outer sample showing smaller Mgb5177 values at given 
a. Similarly, the HgF— a relation clearly differs between the 
two radial samples, with the outer sample having larger HgF 
on average. In this case there is a suggestion of curvature in 
the relation for the core region, with the radial effect more 
pronounced in the low-mass galaxies (log a < 2.1). Finally, 
the Fe5270— a relation shows no discernible change with ra- 
dius. 

Although in principle Figure|2]could form the basis of an 
radial-trends analysis, this method would be quite sensitive 
to the boundaries adopted for the radial bins. Instead, we 
prefer a more flexible approach, treating R c \ as a continuous 
variable, and including an extra term in the index-cr rela- 
tions, to describe the radial behaviour of the line-strengths. 



3.2 Method 

We aim to describe the measured line-strength indices I, 
through a linear model of the form 



I = ao + a a log a + a z lo; 



-Dang(cZCMB) 
-Dang(c«o) 



a R 



Rcl 
R20C 



(2) 



The first two coefficients, do and a a , represent the index-cr 
slope and zero-point as usual. We assume throughout this 
paper that the velocity dispersion, a, can be regarded as a 
proxy for galaxy mass, and moreover that the relationship 
between a and mass is universal. The additional model co- 
efficients, a z and an, respectively test for correlations with 
redshift, and with radius from cluster centre. The principal 
objective of this paper is to constrain and interpret the aR, 
which encode any systematic variation of the line-strengths 
with location in the cluster, after controlling for the over- 
all trend with mass. Thus, a non-zero value of a R would 
not result from the largest galaxies being in the central part 
of the cluster. Rather, aa 7^ indicates a systematic en- 
vironmental dependence in the properties of galaxies at a 
given mass. In Section 14.51 we discuss the effect of using a 
term proportional to logii^i/ifeoo, instead of the linear form 
above. 

Figure |3] shows the sample extent in each cluster, char- 
acterized by the radius R$o% which encloses 90% of the sam- 
ple galaxies, in terms of physical units (upper panel) and af- 
ter scaling by the virial radius (lower panel) . Given that the 
samples were collected with spectrographs of a fixed field of 
view (diameters 60arcmin at WIYN, 38arcmin at CTIO), 



The NO AO Fundamental Plane Survey - III. 5 




Figure 1. The index-cr data for three illustrative indices, Mgb5177, HgF and Fe5270. The index measurements have been adjusted for 
R c l and cz following the best fitting models from Section although in practice this is a small correction. The galaxies used in the 
final fit are denoted with black points and grey error bars. Grey symbols without error bars indicate objects not used in the fit: asterisks 
are galaxies rejected due to nebular emission lines; triangles are objects rejected because they have cD-like J-band luminosities; crosses 
indicate galaxies removed by the iterative outlier rejection scheme. Note that the fitted objects are plotted last, thus only the outlying 
grey points are visible. The best fitting relation is shown by the thick grey line. 




Figure 2. Running quantiles of the Mgb5I77— cr, HgF— a and Fe5270— a relations, for two extreme ranges of cluster-centric radius. The 
'core' region (upper panels) is defined by R c \/R200 < 0-2, while the 'outer' region (lower panels) has 0.5 < R c \/R200 < 1.0. In each panel 
the pale and dark shaded regions enclose respectively two-thirds and one-third of the data in a sliding sample of 300 galaxies (hence the 
reduced range in a relative to Figure^. The global linear index— a fit is reproduced as the bold line for comparison. There is a clear 
offset of the Mgb5177— a relation between the core and outer samples. For HgF, the shift appears more pronounced at low masses, while 
in Fe5270 little difference is observed. 
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the physical area surveyed in a cluster is generally larger for 
more distant clusters than for those nearby. However, when 
scaled in terms of the virial radius, the sample extent is less 
tightly correlated with redshift. This arises partly because 
the cluster sample includes generally more massive systems 
(larger -R200) at greater distances, and partly because there 
is substantial range in the virial radius -R200 at any given 
distance. The cluster with Rgo%/ -R200 > 3 is A3528B, which 
formally has velocity dispersion 50 kms -1 , based on only six 
galaxies. Although such cases clearly yield unreliable -R200 
values, by definition they contribute very few galaxies to the 
fits. Moreover, we restrict the model to radii R c i/R,200 < 1, 
so that such points do not carry undue influence (this cut 
excludes ~4% of the original sample). 
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The parametrization of the redshift term is chosen to 
match the traditional form of the 'aperture correction' (e.g. 
j0rgensen et al. 1995). Here we adopt 10000 kms -1 for the 
normalizing redshift czo; the redshift cz is the CMB-frame 
velocity of the cluster to which the galaxy is assigned, as tab- 
ulated by Smith et al. (2004). We use the ratio of angular- 
diameter distance D ang (cz), for a fi m = 0.3, Qa = 0.7 cos- 
mology, rather than the redshift ratio, but this makes little 
difference in practice. The a z coefficient allows for stellar- 
population gradients internal to the galaxies (e.g. Mehlert 
et al. 2003) , which translate into redshift correlations as the 
observed aperture in kpc scales with distance. Including the 
redshift term allows the results properly to incorporate un- 
certainties in the aperture corrections, and any other sys- 
tematic effects which vary slowly with redshift. In particu- 
lar, any low-order errors in flux calibration will cause weak 
redshift dependence in the indices with 'wide' definitions, 
e.g. CN1, Fe4668). For the most distant clusters, we not 
only observe a larger physical radius in each galaxy, but 
also we typically sample the cluster population out to larger 
radii. This introduces a covariance between an and a z which 
cannot be ignored if meaningful confidence limits are to be 
determined. 

To summarize the fitting method, our analysis uses only 
the emission-free galaxies within one virial radius, excluding 
those with cD-like J-band luminosities. An iterative outlier- 
rejection scheme is employed, clipping galaxies deviating 
from the current model by four times the total measured 
scatter. The number of galaxies in the fit depends on the par- 
ticular index, since no measurements were made for indices 
contaminated by night-sky lines. In practice, ~3000 galax- 
ies are fit for the blue indices up to and including Mgb5177, 
while < 2000 galaxies are fit for Fe5270 and redder indices. 

Our modelling strategy minimises the residuals in the 
line-strengths, implicitly assuming that the index errors 
dominate over the uncertainty in log a. Such a fit will be 
unaffected by explicit selection in a, or in the other pre- 
dictor variables cz and R c \. In practice, however the galaxy 
sample was initially selected on R-band total apparent mag- 
nitude, and the spectroscopic success rate is higher for galax- 
ies of higher central surface brightness. This could introduce 
a slight bias towards observing younger (hence brighter) 
galaxies at any given mass, but one which should not be 
manifest at any particular radius, and so should not affect 
the principal results of this paper. 
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Figure 3. Radial extent of the galaxy sample in each cluster, 
as a function of redshift. The upper panel shows that the more 
distant clusters in NFPS are sampled to a larger physical radius. 
In the lower panel we show the sample extent after scaling to the 
virial radius of each cluster. 



3.3 Results 

The linear fit results are summarized in terms of the 
index— cr relation coefficients, in Tabled 

For the three indices Mgb5177, HgF and Fe5270, the 
global index— a relations were shown in Figure Q Although 
at first sight the slope of the fitted relations appears by eye 
to be somewhat shallow compared to the data, a careful ex- 
amination confirms that the fit approximately bisects the 
vertical distribution at any given a, and thus is a good pre- 
dictor of the line-strengths, as required. 

Figure [I] presents the results on radial trends for all 
twelve indices. Each panel shows the 'partial residuals' of 
the corresponding index— a relation, i.e. 



51 = 1 



ao + a a log a + a z log 



-Dang(cZCMB) 



-Dang(cZo) 



(3) 



averaged in bins of _R c i. The aim is to isolate the variation 
associated with cluster-centric radius trends after correcting 
for the mass and redshift terms. Comparing to Equation [5] 
it is clear that 51 = o,rR c \, if the linear model holds. The 
figure shows our 2a limits on cir as a shaded wedge, for 
comparison to the binned residuals. If the radial behaviour 
is more complex than a linear trend, this will be revealed 
in a poor match between the partial residuals and the fitted 
an. 
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Figure 4. Index-residual versus radius relations. Solid points with error bars show the weighted mean of the partial residuals 81, defined 
as in Equation |21 binned in annuli of constant width 0.1i?200- The SI represent the difference between the measured index and that 
predicted from the velocity dispersion and redshift only, as a function of cluster-centric radius. The residuals have been scaled by the total 
observed scatter in the index— <r relation, the scale factor being noted in each panel. The open symbol at ij c i=0 refers to centrally-located 
cD candidates, which were not used in the fits. The heavy solid line shows our linear fits, with the 2a range of an indicated by a shaded 
wedge. Finally, the heavy dashed line shows the slope predicted under our best-fitting model of Section 14.31 based on a linear radial 
variation in galaxy age, metallicity and [ct/Fe]. 



8 Russell J. Smith et al. 



Table 1. Coefficients of the single-index fits. The total number of galaxies for which the index is measured is given as N && \. After applying 
selection criteria discussed in the text, the sample for each index is reduced to N ae \. In the course of the fitting, a 4<r residual clipping is 
applied, so that the final fit is for N^ t objects. The coefficients ao,a a ,a z ,an are as in Equation 121 For the coefficient of cluster-centric 
radius, an, the significant (2<r) cases are highlighted with an asterisk. The final column gives the ratio of mass-trend to radius-trend, as 
discussed in the text. 



Index 




iV se l 




"0 


a a 


a z 


an 


ratio aji/a a 


HdF 


3765 


3177 


3146 


+3.362 ±0.135 


-1.352 ± 0.060 


+0.107 + 0.061 


+0.036 ± 0.040 


-0.027 ± 0.029 


HgF 


3765 


3177 


3152 


+2.987 ±0.120 


-2.010 ± 0.053 


+0.172 ± 0.051 


+0.136 ± 0.034* 


-0.068 + 0.017 


Hbeta 


3765 


3177 


3169 


+4.397 ± 0.073 


-1.218 ± 0.032 


+0.260 ± 0.030 


+0.073 ± 0.021* 


-0.060 + 0.017 


CN1 


3765 


3177 


3145 


-0.333 + 0.008 


+0.179 ± 0.003 


-0.024 ± 0.004 


-0.006 ± 0.002* 


-0.033 + 0.013 


Mgb5177 


3716 


3142 


3139 


-1.849 + 0.085 


+2.780 + 0.037 


-0.227 + 0.036 


-0.203 ± 0.025* 


-0.073 ± 0.009 


Fe4668 


3765 


3177 


3165 


-4.090 ± 0.248 


+4.735 ± 0.110 


-1.381 ± 0.112 


-0.163 ± 0.073* 


-0.035 + 0.016 


Fe4531 


3765 


3177 


3160 


+1.537 + 0.110 


+0.779 ± 0.049 


-0.260 ± 0.047 


+0.072 ± 0.032* 


+0.092 ± 0.042 


Fe5270 


2152 


1824 


1822 


+1.797 + 0.087 


+0.468 ± 0.038 


+0.077 ± 0.042 


+0.034 ± 0.028 


+0.072 + 0.061 


Fe5335 


2291 


1966 


1962 


+1.368 + 0.096 


+0.554 ± 0.043 


-0.053 ± 0.037 


+0.024 ± 0.031 


+0.044 ± 0.056 


Fe5406 


3071 


2583 


2580 


+0.990 ± 0.063 


+0.309 ± 0.028 


-0.102 ± 0.023 


+0.017 + 0.018 


+0.056 ± 0.058 


Fe5709 


1831 


1524 


1522 


+1.133 + 0.055 


-0.143 ± 0.024 


+0.007 ± 0.022 


+0.039 ± 0.017* 


-0.270 + 0.124 


Fe5782 


1047 


888 


887 


+0.211 + 0.065 


+0.239 ± 0.029 


+0.091 ± 0.047 


+0.045 ± 0.022* 


+0.188 + 0.094 



For our three illustrative indices, we show the error el- 
lipses in Figure demonstrating the modest covariance be- 
tween an and the other coefficients. The correlation coeffi- 
cient for the an and a z parameters is —0.32 < r < —0.24 
with larger values for the red Fe lines which are available 
over a narrower range of redshift due to sky-line contami- 
nation and the more limited spectral extent at CTIO. The 
cause of this effect is the correlation of sample radial extent 
with redshift, as discussed above. There is also a small co- 
variance between an and a a (0.13 < r < 0.19), in the sense 
that fitting a flatter slope with respect to mass causes a 
steeper radial gradient. This cross-talk arises because sam- 
ple galaxies near the cluster centre have a larger a on average 
than those further out (~0.1dex out to the virial radius), 
which could be due to true mass segregation in the cluster, 
but could also be an artefact of greater target density in 
the centre, coupled with our prioritizing brighter targets for 
spectroscopic follow-up. In practice, the level of covariance 
is small enough that reasonable changes (~10%) in a a do 
not significantly influence the radial trend. 

Taken jointly, the overall significance of the _R c i trends 
is very high. With 12 indices we would expect of order four 
of the an coefficients to be non-zero at la by chance, includ- 
ing perhaps one at 2a. In fact, the Mgb5177 trend is signif- 
icant at ~8er, while several others are at > 3<r. Relative to 
a model in which the indices are uncorrelated with cluster- 
centric radius, we formally have \ 2 ~ 126 on 12 degrees of 
freedom. We can therefore already exclude the null hypoth- 
esis and conclude with high degree of confidence that the 
line-strengths of red-sequence galaxies in the NFPS dataset 
exhibit signicant variations as a function of cluster location. 

It is interesting to consider ratio the an/a a , which com- 
pares the effect on a given index of increasing the galaxy 
mass with the effect of increasing the cluster-centric radius. 
A simple hypothesis would be that galaxies at larger radius 
have the the same stellar population properties as more cen- 
tral galaxies of different mass. In this case, the stellar pop- 
ulations would form a one-parameter family; for example, 
this could arise from a single evolutionary process which is 
switched off at a different redshift, dependent on both a 
and Rci- A one-parameter model of this form would pre- 



dict that all indices have the same an/a a ratio within the 
errors. In fact, as shown in Tableland Figure |S| there is 
clearly a variation in this ratio among the twelve fitted in- 
dices. Specifically, in the Balmer indices (HdF, HgF, Hbeta) 
and the indices Mgb5177, CN1 and Fe4668, the effect of in- 
creasing the radius is equivalent to that of decreasing the 
mass. However, for most of the Fe indices, the effect of in- 
creasing the radius is similar to increasing the mass. We 
can infer from this that the processes driving the radial sys- 
tematic behaviour are distinct from those which imprinted 
the index — a relations themselves 3 . This is quite a general 
conclusion, and does not require any model-dependent inter- 
pretation of the indices. To anticipate the results of the next 
section however, note that the indices with an/a a < are 
those associated with age and [a/Fe] in the TMBK models, 
while the indices with an/a a > trace mainly the over- 
all metallicity [Z/H], Thus, we may expect that the scaling 
of the three parameters with cluster-centric radius will dif- 
fer qualitatively from the mass-scaling relations. A similar 
conclusion can be reached from Figure which shows two 
classical 'two-index' diagrams. In particular, the diagram of 
Mgb5177 vs (Fe) shows that the mass-trend followed by the 
outer galaxies is offset from that of the inner objects. The 
effect of changing cluster-centric radius is not parallel to the 
effect of changing the galaxy mass, again suggesting that at 
least two distinct mechanisms are at work in determining 
the average stellar population at a given mass and radius. 



4 INTERPRETATION OF THE GRADIENTS 

In this section, we analyse the single-index fit results in 
terms of the TMBK spectral synthesis models. The models 
predict the values of each index, over a grid in the three- 
parameter space of age, total metallicity ([Z/H]) and a- 
element to iron-peak element abundance ratio ([a/Fe]). Fol- 
lowing the differential method introduced by Nelan et al. 

3 Moreover, the range in an/a a argues against an explanation in 
terms of an environmentally-driven systematic shift in the veloc- 
ity dispersion, e.g. a non-universal mass-cr relation. 
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Figure 5. Error ellipses for the single-index fits, for the example cases of Mgb5177, HgF and Fe5270. The upper row shows the 1,2,3-ct 
confidence regions for the coefficients of cluster-centric radius and redshift in Equation [5] In the lower rows, the equivalent regions are 
shown for the coefficients of radius and log a. Dashed lines indicate zero for reference. 
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Figure 6. Ratio of radial-slope to index— a slope, i.e. a^/a-a- The shaded boxes indicate the lex limits on this ratio as given in Table HI 
The vertical bars divide the indices into subsets broadly reflecting their responses to age, metallicity and a/Fe ratio, in the models of 
TMBK. The significant index-to-index variation in aji/a a argues for an environmental mechanism which is distinct from the overall mass 
trend (see text). 
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Figure 7. Classic two-index diagrams. The left panel shows Hbeta versus [MgFe'] (defined as in Thomas et al. 2003), while the right 
panel shows Mgb5f77 versus (Fe> (=±(Fe5270+Fe5335)). In each case, the five pairs of boxes indicate five bins in velocity dispersion; 
these bins are defined by a ranges as in Table 8 of Nelan et al., the mean a being indicated in the plot. For each bin, the dark-shaded 
box indicates the la error on the weighted mean line-strengths for galaxies with R c \/R200 < 0.3 (approximately the median R c i), while 
the light-shaded box shows the same for 0.3 < R c \/R200 < 1-0. Bold vectors show the effect of f5% increases in age, metallicity and 
[a/Fe], on each pair of indices, from the models of TMBK. 



(2005), we attempt to reproduce the observed index- versus- 
radius trends (i.e. the values of or) in terms of underlying 
radial trends of the three model parameters. 



4.1 Parameter responses in the TMBK models 

In principle, given line-strength measurements of sufficiently 
high signal-to-noise, and assuming the calibration of the 
models and data to be secure, it should be possible to 
convert the index measurements for a single galaxy into 
determinations of age, [Z/H] and [a/Fe] for that object. 
This 'grid-inversion' approach has traditionally utilised a 
small number of indices such as Mgb5177, Hbeta and 
Fe5270/Fe5335 (e.g. Trager et al. 2000; Kuntschner et al. 
2001, 2002; Poggianti et al. 2001; Thomas et al. 2005), but 
can be readily extended to a larger index set (Proctor & 
Samson 2002; Proctor et al. 2004). 

In practice however, the observed line-strengths are sub- 
ject to large random errors, so that in general some measure- 
ments scatter outside of the model grid, and because the 
model grids are tilted with respect to the observed quanti- 
ties, errors in the derived population parameters are corre- 
lated, complicating their interpretation (Kuntschner et al. 
2001; Thomas et al. 2005). Beyond the problems associated 
with random errors, which can be addressed by considering 
composite spectra (Bernardi et al. 2006) or averaged indices 
(Nelan et al. 2005), there are also systematic calibration un- 
certainties in both the data and the models. This is often 



expressed by the observation that, for instance, the overall 
position of the classical Mgb5177 versus Hbeta diagnostic 
grid is less secure than the relative differences in predictions. 
These considerations motivated us to introduce a differen- 
tial analysis of the index-<r relations, (Nelan et al. 2005). 
Here, we recap this approach, and then extend the method 
to describe the cluster-centric trends in the same manner. 

We use the TMBK models to predict the linear 're- 
sponses' of each index to each of the driving parameters, 
while holding the others constant. Thus for Mgb5177, we 
determine the following responses: 



n,Mgb5177 
-'''age 



R 



d Mgb5177 
d log Age 



M g b5i77 _ d Mgb5177 
d [Z/H] 



Z/H 



and 



R 



Mgb5177 
a/Fe 



d Mgb5177 

d [a/Fe] 



Such a description is valid where the model grids are ap- 
proximately parallel and of constant gradient. This is the 
chief simplifying assumption, and also the main limitation of 
our method. The responses are computed relative to a fidu- 
cial population, and over a subset of the TMBK parameter 
space, roughly corresponding to old, metal-rich, a-enhanced 
ellipticals (Table |2J • For example, we determine R^^ 5177 
over the range 0.00 — 0.35 in [Z/H], holding the age at 8Gyr 
and [a/Fe] at 0.2. (Strictly, as in Nelan et al., the responses 
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Table 2. Fiducial model and the parameter ranges used to com- 
pute parameter responses from the TMBK stellar population 
models. 



index 
log age 



plo gagC+ ~ 



index r>[Z/H] 
'Z/H] cr 



i/Fc] 



(4) 



Parameter 


fiducial 


rang< 




Age (Gyr) 


8.0 


3.0 - 


15.0 


[Z/H] 


0.00,0.35 


-0.33 - 


0.67 


[a/Fe] 


0.20 


0.00 - 


0.50 



Table 3. Index responses to population parameters. i? par is the 
change in the index per decade in parameter 'par', estimated over 
the range given in Table [3] and holding the other parameters 
at their fiducial values. As an estimate of the errors in the re- 
sponses, Spar is the rms scatter in i? par as computed over pairs 
of models within the parameter range. The quantities arc for the 
line-strengths expressed in Angstroms, except for CN1 in magni- 
tudes. The divisions in the table reflect the relative responses, as 
in Figure 151 



Index 


Rage 




Ret /Fc 


•Sage 


•Sz/H 


Sct/Fe 


HdF 


-1.921 


-1.393 


1.889 


0.445 


0.414 


0.752 


HgF 


-2.803 


-2.223 


1.091 


0.599 


0.391 


0.214 


Hbeta 


-1.134 


-0.603 


0.266 


0.211 


0.189 


0.053 


CN1 


0.098 


0.177 


0.052 


0.046 


0.032 


0.019 


Mgb5177 


1.547 


2.914 


2.026 


0.349 


0.498 


0.634 


Fe4668 


1.748 


7.720 


0.608 


0.918 


0.740 


0.707 


Fe4531 


0.812 


1.565 


-0.940 


0.243 


0.223 


0.178 


Fe5270 


0.660 


1.488 


-1.305 


0.242 


0.175 


0.230 


Fe5335 


0.577 


1.729 


-2.078 


0.179 


0.461 


0.598 


Fe5406 


0.437 


1.181 


-1.157 


0.145 


0.208 


0.237 


Fe5709 


0.082 


0.573 


-0.312 


0.048 


0.033 


0.019 


Fe5782 


0.104 


0.527 


-0.344 


0.064 


0.036 


0.027 



for age and [a/Fe] are an average of those for fiducial metal- 
licities [Z/H]=0 and those for [Z/H]=0.35.) 

The complete set of response vectors, as used for the 
analyses described below, is summarized in Table |3] The ta- 
ble also indicates the scatter, S, in the response, computed 
over all pairs of TMBK models within the range specified 
in Table |2] This quantity serves as a conservative estimate 
of the error introduced by forcing a linear description of the 
grids. In Figure |§1 we show the age and [a/Fe] responses 
of the twelve indices relative to their [Z/H] response. This 
figure neatly separates the three classes of index described 
in Section |21 and in particular shows why we classify Fe4668 
with Mgb5177 and CN1 rather than with the other Fe in- 
dices. It also shows why the behaviour of HdF is expected 
to differ from that of Hbeta and HgF, as discussed later. 

As in Nelan et al., we model the observed index— a 
slopes in terms of underlying linear trends of log(age), [Z/H] 
and [a/Fe], as a function of log a, using the response vectors 
to translate between the observed and underlying relation- 
ships: 



(la 



dl 
dlogCT 

dl d log age 



d log age d log a 



dl g[Z/Hj dl 
+ <9[Z/H] ' d log a + d[a/Fe 



d[a/Fe] 
d log a 



where the change of each parameter with log a has been 
written as P£ al for brevity. 

We can extend the Nelan et al. method to the radial 
trends, by writing an equivalent expression to describe the 
an coefficients in terms of underlying changes of age, metal- 
licity and a-abundance as a function of cluster-centric ra- 
dius, R ci : 



log age . rjindex D [Z/H] jyindex t-jL ( 
+ -fMZ/Hl -rj?_, + n[ a /p c ]-r R : 



:/Fe] 



(•>) 



With at least three indices having different response 
signatures relative to the three population parameters, we 
can attempt to invert these equations for the elements of 
Pa and Pr^ , and hence determine how age, metallicity and 
a-abundance scale with mass and with location in the clus- 
ter. A minimal set of indices for this purpose might be HgF, 
Mgb5177 and Fe5270. The first of these provides age sensi- 
tivity while the remaining two trace metallicity, with differ- 
ing dependency on [a/Fe]. In what follows, we use not only 
these three but all of the twelve indices chosen in Section |2| 
to over-constrain the solution. We cast the fit in terms of 
predicting the coefficients an, with the free parameters 
to be determined. In the fitting stage, we assume the R^ ex 
are known precisely, and weight the data only according to 
the errors on the measured clr. 

In summary, we emphasize that our strategy for deter- 
mining the parameter gradients P is explicitly a two-step 
process: first, we fit the line-strength data for the -~3000 
galaxies to determine the radial coefficients of the indices; 
second, we model these twelve coefficients an by fitting to 
the model responses, reducing the information to just three 
parameters. If a good fit can be obtained, the slopes of age, 
[Z/H] and [a/Fe] with radius provide an economical descrip- 
tion of the observed trends. A weakness of this approach is 
that it implicitly assumes the errors on the an are not corre- 
lated with each other, such that each index slope constitutes 
an independent constraint on the model. A more sophisti- 
cated approach would attempt to combine the two steps into 
a single fit. The present data probably do not justify such a 
treatment, and we prefer to retain the intermediate step for 
clarity of presentation. 



4.2 Stellar population trends along the a sequence 

First, we revisit the mass scaling relations, as earlier dis- 
cussed by Nelan et al. (2005). Figure [5] shows the confidence 
intervals on the quantities P^ ar , derived from the index— a 
slopes a a , as in Table We derive the following scalings: 



A U 

Age oc a 



Z/Hoc a 



a/Fe oc g 



(6) 



confirming significant trends of all three parameters with 
a, so that on increasing the velocity dispersion, the stel- 
lar populations become older, more metal rich, and more 
a-enhanced. As indicated in the figure, there is a strong de- 
generacy between the age and metallicity scalings, such that 
the errors on these parameters are correlated. Comparison 
with the trends reported by Nelan et al. shows that despite 
several differences in the sample construction and analysis 
methods (notably, the presence of additional terms in the 
index— a model), the a a coefficients yield scaling relations 
consistent with the previous results, at the ~la level. 
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Figure 8. Linear responses of indices in the TMBK models, over 
the parameter range appropriate to early- type galaxies. The age 
and a-abundance response are each normalized by the metallicity 
response, to reduce the information to a two-dimensional format. 
The quadrants, formed by arbitrarily-placed vertical and hori- 
zontal lines, divide the indices according to the parameters which 
they best constrain, as in Figure 151 



Figure 1101 provides a visual impression of the success 
of Equation |S|in reproducing the measured index— a slopes 
(cf. Fig. 10a of Nelan et al.). We plot the measured slopes 
for each index on the horizontal axis, normalized by the 
measurement error. Thus, the position of each index on the 
horizontal axis reflects its significance, and the horizontal 
errors are unity by definition. The vertical axis gives the 
predicted index-er slopes, under the scalings of Equation [(J 
also normalized by the error in the measurement. The er- 
rors in the predictions are derived from the scatter in the 
responses over the relevant subregion of the model grid, i.e. 
the Spar given in Table |3f. For an adequate solution, the 
error boxes for all indices would lie along the line of unit 
slope. Outliers from this line highlight indices with a a co- 
efficients not well described by the model. For instance, as 
discussed by Nelan et al., Fe5709 decreases with increasing 
cr, although our model would predict it to increase, driven 
primarily by [Z/H]. 



4.3 Cluster-centric trends in stellar populations 

We now turn to the main objective of this paper, which is to 
determine which aspects of the stellar populations are dif- 



4 For the mass-scalings discussed here, errors in the predictions 
dominate over the errors in the observed slopes. However, the 
prediction errors are only indicative, and are not independent 
from index to index. Thus we do not attempt to derive a formal 
goodness-of-fit in this case. In contrast, for the radial trends anal- 
ysed later, the observed errors are comparable to the prediction 
errors, and a \ 2 can De quoted with greater justification. 



ferent in the cores of clusters compared to the outer regions, 
and are responsible for driving the radial dependencies of the 
index — o residuals. Unlike many previous works (Guzman et 
al. 1992; Carter et al. 2002), our approach explicitly allows 
for simultaneous variations in age, total metallicity and a- 
abundance ratio. 

Figure [TTI shows the confidence intervals on the quan- 
tities P^ a , derived from the coefficients an. The key result 
is that both age and [a/Fe] are found to depend on cluster- 
dentric distance, while no there is no significant trend of 
[Z/H]: 

d log Age 



dR cl 

d[Z/H] 
di? c i 

dJo/Fe] 
di?d 



-0.072 ±0.017 
+0.007 ±0.010 
±0.048 ±0.009 



(7) 



Thus, galaxies at 1 -R200 are on average 15% younger and 
10% less a-enriched than galaxies of the same mass located 
in the core of the cluster. Considering the one-parameter 
confidence intervals, the significance of the age and [a/Fe] 
variations are above the 4<r and 5cr levels, respectively. 

The success of this model, in reproducing the observed 
cluster-centric trends in the indices, can be judged from Fig- 
ure^] which compares the measured and predicted an for 
each index. The figure confirms that the model indeed ac- 
counts for the decrease of Mgb5177, CN1 and Fe4668, and 
also the increase of Hbeta and HgF with increasing radius. 
The case of the remaining Balmer index HdF is instruc- 
tive. The age-sensitivity of HdF would lead to a positive 
slope with radius, if the age-radius trend acted alone. How- 
ever, compared to Hbeta and HgF, the HdF index also re- 
sponds strongly to [a/Fe], as seen in Figure |H] At large 
radius, the HdF is strengthened by the younger ages, but 
weakened by the lower [a/Fe]. Our model predicts the ef- 
fects of age and [a/Fe] gradients will compensate each other 
in the case of HdF, successfully reproducing the observed 
fairly flat behaviour. The Fe-sensitive indices are predicted 
to have weak trends, as observed. The model predictions for 
each index are also over-plotted on the partial residual dia- 
grams in Figure 2] for comparison to the linear fits for an, 
and to the binned data. It is seen that only in the case of 
Fe4531 is the model prediction outside of the 2a observa- 
tional limits. In contrast to the equivalent diagram for the 
mass-scalings (Figure ED, the measurement errors are com- 
parable to the errors in the predictions. Taking both into ac- 
count, the model has x 2 = 10-6 for nine degrees of freedom 
(twelve parameters less three fitted coefficients), and thus is 
formally an acceptable fit. Hereafter, we refer to the solution 
given in Equation [7| as the 'default model'; its parameters 
are summarized in Line 1 of Table 0]for reference. For com- 
parison, Line 2 of the table gives the \ 2 = 126 obtained 
for the null model (no dependence of stellar populations on 
cluster-centric radius). 

To illustrate further the influence of particular indices 
in our preferred solution, we consider whether the an can be 
adequately fit by 'restricted' models, in which only one or 
two of the stellar population parameters are allowed to vary 
with R c \ . The results are shown in Figure 1131 where each 
panel is equivalent to Figure E3 but now for the best-fitting 
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Figure 9. Constraints on variation of TMBK population parameters with velocity dispersion. The ellipse boundaries indicate the 68%, 
95% and 99% confidence levels in the two parameters taken jointly. The black crosses show the results derived by Nelan et al. (2005). 



Table 4. Cluster-centric trends of age, [Z/H] and [a/Fe], The various lines record the results of tests discussed in the text, as follows: 
Line 1 gives the results from the default three-parameter model (Section l4.3l : Line 2 gives the x 2 for the null hypothesis for comparison; 
Lines 3-5 show the best-fitting two-parameter models; Lines 6—8 show the best-fitting one-parameter models; Lines 9-10 investigate 
the influence of morphology on the solution (Section 14.41 ; Lines 11-21 test for systematic errors by fitting subset of the galaxy sample 
(Section l4~5t : 







log(Age) gradient 


[Z/H] gradient 


[a/Fe] gradient 


X 2 /dof 


notes 


1 


Default 


-0.072 ±0.017 


±0.007 ±0.010 


-0.048 ± 0.009 


10.6 / 9 


default solution (Equation 171 


2 


No gradients 








126.1 / 12 


null model 


3 


age & [a/Fe] 


-0.063 ± 0.009 




-0.048 ± 0.008 


12.0 / 10 


no radial [Z/H] trend 


1 


[Z/H] & [a/Fe] 




-0.031 ±0.010 


-0.043 ± 0.014 


40.1 / 10 


no radial age trend 


5 


age & [Z/H] 


-0.059 ± 0.033 


±0.004 ± 0.021 




50.0 / 10 


no radial [a/Fe] trend 


6 


Age only 


-0.054 ± 0.017 






53.1 / 11 


no [Z/H] or [a/Fe] trends 


7 


[a/Fe] only 






-0.037 ±0.019 


75.5 / 11 


no age or [Z/H] trends 


8 


[Z/H] only 




-0.027 ±0.012 




79.7 / 11 


no age or [a/Fe] trends 


9 


B/T > 0.5 


-0.040 ± 0.023 


-0.013 ±0.014 


-0.048 ±0.012 


8.9 / 9 


bulge-dominated 


10 


fit B/T 


-0.066 ± 0.024 


±0.002 ± 0.015 


-0.046 ± 0.013 


12.5 / 9 


control for B/T in index fits 


11 


logo- > 2.10 


-0.051 ±0.013 


-0.001 ±0.008 


-0.054 ± 0.007 


5.9 / 9 


excludes ~ -| of galaxies 


12 


logo- < 2.24 


-0.096 ±0.019 


±0.004 ± 0.012 


-0.048 ± 0.010 


7.1 / 9 


f 

excludes ~ ^ of galaxies 


13 


DX clusters 


-0.100 ± 0.021 


±0.018 ±0.013 


-0.051 ± 0.011 


9.7 / 9 


coincident cD and X-ray peak 


14 


cz < 13800 kms" 1 


-0.061 ± 0.029 


±0.014 ±0.018 


-0.045 ± 0.014 


14.1 / 9 


clusters within median depth 


15 


cz > 13800 kms" 1 


-0.078 ± 0.033 


-0.001 ±0.021 


-0.047 ±0.018 


22.8 / 8 


clusters beyond median depth 


16 


Dec. > 0° 


-0.096 ± 0.026 


±0.015 ±0.015 


-0.035 ± 0.013 


10.3 / 9 


mostly WIYN 


17 


Dec. < 0° 


-0.062 ± 0.014 


±0.009 ± 0.009 


-0.068 ± 0.008 


4.5 / 9 


mostly CTIO 


IS 


No Mj cut 


-0.081 ± 0.015 


±0.011 ±0.009 


-0.051 ± 0.008 


9.2 / 9 


includes cDs 


19 


No emission cut 


-0.073 ± 0.017 


±0.003 ± 0.010 


-0.051 ± 0.008 


9.5 / 9 


No cut in [O in] or H 


20 


R cl < 2.0i?200 


-0.068 ±0.017 


±0.007 ±0.010 


-0.033 ± 0.009 


16.7 / 9 


increased radial range of fit 


21 


R cl < O.5R200 


-0.139 ± 0.030 


±0.028 ±0.019 


-0.072 ± 0.015 


13.1 / 9 


reduced radial range of fit 



restricted model of the given type. It is important to note 
that this is not a projection of the default fit into a reduced 
parameter-space; i.e. the [Z/H]-vs-radius slope of the model 
considered in Figure I13H (for instance) is not the same as 
that in Figure 1121 As expected from the coefficients of the 
default model, removing the radial dependence of [Z/H] does 
not significantly compromise the fit (Figure [F3"k and Line 3 
of Table QJ. By contrast, we find a very poor fit for all of 
the models which exclude either the age or [a/Fe] variations. 
In particular, the best-fitting model without cluster-centric 
age variations (Figure 112b and Line 4) cannot reproduce 



the increase of Hbeta and HgF with radius. Models with- 
out cluster-centric [a/Fe] variations (Figure lTSl : and Line 5) 
cannot account for the strong negative trend of Mgb5177, 
and also predict a weak negative trend for the Fe indices, 
rather than the positive trends observed. Finally, we obtain 
extremely poor fits for models which attempt to account 
for the radial trends in terms of a single driving parameter 
(Figure lF3l i.e.f and Lines 6-8). Taking these results at face 
value, any single-parameter model, e.g. a pure age trend (as 
in Guzman et al. 1992) or a pure metallicity trend (as in 
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Figure 10. Predicted versus observed indcx-cr slopes in the best 
fitting model for the trends with mass. Both the model predictions 
and the observed slopes are scaled by the error in the observed 
slope. The shaded box shows the lcr errors, where the errors in the 
predictions are derived from the errors in the responses fTablcl^l. 
Indices which deviate significantly from the unit-slope line have 
slopes poorly reproduced by the model. 

Carter et al. 2002), is bound to yield inconsistent results 
when applied to multiple indices. 

4.4 Influence of morphological trends 

A very pertinent concern is whether the observed gradients 
truly represent changes in stellar populations of otherwise 
similar galaxies, or whether instead they reflect a radial vari- 
ation of the morphological mixture in clusters. Recall that 
the galaxy sample was selected using colour and magnitude 
information only, without reference to galaxy morphology 
(Smith et al. 2004). Additionally, for the present analysis, 
we imposed a selection against objects with nebular emission 
lines in the (central) spectra. Thus although star-forming 
discs should be excluded, is is very likely that our sample 
includes some bulges of early- type spirals. If such systems 
have different stellar populations from pure ellipticals of the 
same mass, and if the fraction of spirals in the sample in- 
creases in the cluster outskirts, then this might be manifest 
as radial trends of the kind observed. 

At present, quantitative morphological information is 
available for approximately half of the NFPS galaxies, and 
of these around 70% have bulge-dominated morphologies 
(bulge-to-total light ratios, B/T, greater than 0.5). One sim- 
ple test for morphology-driven trends is to fit only these 
bulge-dominated galaxies. In this case, the detection of an 
age trend is reduced to a confidence level below 2cr, although 
the [a/Fe] trend remains significant (Line 9 of Tabled . How- 
ever, the errors are large after this sub-sampling, and the de- 
rived gradients are also less than 2<r away from our default 
solution. If the difference is real, this might mean that part 



of the radial trends are driven by changing morphological 
composition with radius. However, it could also result from 
a true radial age effect (i.e. at fixed morphology), but due 
to a mechanism that is more effective in galaxies with discs. 

An alternative way to test for morphology effects is 
to allow an additional term proportional to B/T in the 
index— a fits. To first order, this term should absorb any 
systematic dependence on morphology, so that the R c \ term 
describes the radial effect at fixed bulge-fraction as well as 
fixed mass. Fitting this more complex model to the subsam- 
ple with measured B/T indeed reveals a strong morphology 
dependence in the index data. The pattern of effects is dif- 
ferent again from either the mass sequence or the radial 
trends. Specifically, the Balmer lines and Fe-dominated in- 
dices are all decreased for more bulge-dominated systems, 
while the a-dominated indices Mgb5177 and CN1 increase 
with greater bulge-fraction. We can quantify this effect in 
terms of stellar population parameters, by modelling the 
morphology coefficients a B / T in the same way as we treated 
a a and an. We find that for a shift of one unit in B/T (i.e. 
going from pure disc to pure bulge), we increase log(Age) 
by 0.176 ± 0.026, decrease [Z/H] by 0.097 ± 0.016 and in- 
crease [a/Fe] by 0.052 ± 0.013. In summary, galaxies with 
larger bulges are older, less metal rich and more a-enhanced 
than diskier objects of the same velocity dispersion. At face 
value, this is dissimilar from the conclusions of Thomas & 
Davies (2006), who find that early-type spiral bulges have 
similar stellar populations to ellipticals of the same velocity 
dispersion. 

Having controlled for the morphology effects, we can 
now assess their impact on the cluster-centric trends. Per- 
haps surprisingly, the corrected radial trends (Table 2] Line 
10) are little changed relative to our default solution. This 
robustness arises because the morphology-radius relation is 
rather shallow for our sample, which was already selected 
to have red colours and no nuclear emission. From the core 
of the cluster out to the virial radius, the average B/T is 
reduced only by 0.073 ± 0.014. Combining this figure with 
the above morphology effects yields a 'correction' which is 
consistent with the difference between the default and the 
morphology-controlled solutions. In passing we note that the 
mass-scaling trends are also not strongly affected by the in- 
clusion of the B/T term, despite a strong morphology-<r re- 
lation (B/T increasing by 0.39 ± 0.03 per decade in a). In 
this case, when multiplied by the morphological effect on the 
stellar populations above, the net effect is an order of mag- 
nitude smaller than the mass trend at constant morphology. 
For instance, for one decade in a, morphology drives an age 
change of 0.18 x 0.39 = 0.07 in the log, compared to the 
derived relation Ageoc <t~ ' 7 . 

We conclude that although red-sequence galaxies of dif- 
ferent morphologies indeed seem to have systematically dif- 
ferent stellar populations, this has little impact on the re- 
sults derived here for the cluster-centric trends, or for the 
mass-scalings. An investigation of morphology effects on the 
NFPS colour-magnitude relations is presented by Stevenson 
et al. (2006). A full synthesis of morphological and stellar 
population information must await completion of the B/T 
measurements, and will be addressed in a future paper. 
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Figure 11. Constraints on variation of TMBK population parameters with cluster-centric radius. The ellipse boundaries indicate the 
68%, 95% and 99% confidence levels in the two parameters taken jointly. 
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Figure 12. Predicted versus observed slopes in the best fitting 
model for the trends with cluster-centric radius. The diagram is 
equivalent to Figure [TOl for the mass-scalings. 

4.5 Further robustness tests and data subsamples 

In this section we continue to address the robustness of the 
results, by testing the fits obtained after sub-sampling the 
input data set. The results are summarized in the lower sec- 
tion of Table H 

First, motivated by the apparent mass-dependence of 
the HgF trends seen in Figure we test whether the re- 
sults are strongly driven by the low-mass objects. Restrict- 
ing the sample to galaxies with fairly high velocity dispersion 
a > 125 kms -1 removes around one third of the objects. Fit- 
ting to the reduced sample yields a slightly reduced age gra- 
dient, but the solution remains marginally consistent with 
our default fit (Line 11 of Table 0. Similarly, restricting to 
a < 175 kms -1 yields a marginally stronger radial age gra- 
dient (Line 12); the [a/Fe] and [Z/H] trends are unaffected 
by cuts in a. Despite sharing 50% of galaxies in the inter- 



mediate mass range, the two subsamples differ from each 
other at ~2<r in the recovered age gradient. Thus there is 
some evidence that the age trends are more pronounced for 
galaxies of smaller mass. Because of the tendency towards 
more disc-dominated objects at low <r, this test is not in- 
dependent of the B /T cut in Section 14.41 Note also that 
because the average stellar populations vary strongly with 
mass, we cannot make arbitrarily severe cuts on a without 
potentially invalidating the response vectors adopted. 

A potential concern is that some clusters have ambigu- 
ous or poorly-determined centres (Smith et al. 2004), such 
that, in the worst cases, the R c \ values used in the fits will 
not be a good proxy for galaxy environment. If we use only 
the galaxies in clusters with the most secure centres ('DX', 
based on a dominant galaxy coincident with an X-ray emis- 
sion peak) , we obtain a marginally stronger and more signif- 
icant age trend (Line 13). This suggests that clusters with 
less well-determined centres are diluting the trends in the de- 
fault solution, either because the radii are less meaningful or 
else because more relaxed clusters have intrinsically stronger 
gradients. In this context, we note that in some well-sampled 
systems such as instance Abell 3558, the radial gradients in 
Mgb5177, etc, can be detected for individual clusters. This 
in itself confirms that the results obtained for the full sample 
are not simply an artefact of stacking the clusters. 

In more distant clusters, we typically sample more lu- 
minous galaxies on average, and to somewhat larger cluster- 
centric radius. Moreover, any instrumental systematic error 
affecting data over a particular wavelength range will have 
an impact on different indices at different redshifts. A simple 
test for such effects is to divide the sample into halves ac- 
cording to host cluster redshift, as in Lines 14-15. The two 
cluster sets yield radial trends which are consistent with 
each other, and with the default solution, within the en- 
larged errors. An alternative division of the cluster set is 
into northern and southern celestial hemisphere, which are 
dominated, respectively, by data obtained at the WIYN and 
CTIO telescopes. Systematic offsets between index measure- 
ments at one telescope relative to another have been investi- 
gated at the data level by Nelan et al. (2005). Any residual 
offsets could in principle lead to radial trends in the data 
since the WIYN field of view was wider than the equivalent 
field of view at CTIO. Fitting the northern and southern 
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cluster sets separately, we find a weaker [a/Fe] gradient in 
the northern data than the southern data (at the ~2cr level; 
Lines 16-17). Each is marginally inconsistent with the de- 
fault solution. The age gradients are consistent between the 
subsamples. 

Next we relax some of the selection criteria used to de- 
fine the galaxy sample, to check how these affect the results. 
Removing the cut on J-band luminosity, allowing the cDs 
(and other large galaxies) to re-enter into the fits, there is 
no significant effect on the results (Line 18). If we allow 
galaxies with nebular emission to contribute to the fits, the 
solution is again almost unaffected (Line 19) . This surprising 
stability arises partly because the iterative clipping process 
efficiently removes many of the galaxies with emission, and 
partly because we are using multiple Balmer indices, with 
the more robust HgF and HdF balancing the more easily 
contaminated Hbeta. 

Finally we test the effect of changing the outer radial 
cut. We noted in Section r3.2l that a few clusters have unreli- 
able -R200 due to poorly-determined cluster velocity disper- 
sions. In such clusters, the R c i/R200, are subject to system- 
atic errors. To reduce the influence of galaxies in these clus- 
ters, we limited the fit to i? c i< R200 • Replacing some of these 
galaxies, and fitting now over a larger range, i? c i< 2i?2oo, 
we derive a consistent age trend, while the trend of [a/Fe] 
is marginally weakened (Line 20). If instead we restrict the 
range still further, to R c i< 0.5-R200, both the age and [a/Fe] 
trends become stronger (by ~2<r and ~1.5<r respectively), al- 
beit with errors 50% larger (Line 21). This behaviour can be 
discerned in the radial residuals diagram (Figure [IJ, where 
a number of the indices show a stronger trend in the in- 
ner few bins, and are flatter further out. Motivated by the 
apparently stronger gradients in the inner regions, we have 
experimented with an alternative parametrization of the ra- 
dial trends, using a term proportional to log(i? c i/i?20o)- This 
requires an ad hoc treatment of the galaxies at very small 
Rd, to prevent them from exercising excessive leverage on 
the fits. If we place the inner galaxies at a minimum radius 
of 0.01 -R200, we recover an age trend of 12±2% per decade 
in R c \/R200 and an a/Fe trend of 7±1% per decade. These 
results are, however, sensitive the treatment adopted for the 
innermost galaxies. 

In conclusion, among all of the data subset tests per- 
formed above, significant trends towards younger ages and 
lower [a/Fe], at larger cluster-centric radius, are recovered 
in every case. In addition to demonstrating the robustness 
of the cluster-centric variations, these tests have hinted at a 
number of possible departures from the assumptions of our 
modelling. First, there is some evidence that the cluster- 
centric trends are steeper in the cluster core than further 
out; alternative parametrizations of the R c \ term could in 
principle yield improved fits. Second, the trends appear to 
be stronger for low-mass and/or disky galaxies than for high- 
er bulge-dominated systems. Our model (Equation |5J does 
not naturally account for such an effect, since it assumes lin- 
ear addition of mass and radius contributions, with no cross 
term (i.e. no term proportional to the product 7? c ilog<r). 
However, the results of the robustness tests might be point- 
ing towards real differences in the impact of the environment 
on different galaxy types and masses. 



5 DISCUSSION 

In this section we first compare our results with previous 
work on trends in red-sequence galaxy properties within 
clusters. We also review the distinct approach of comparing 
cluster galaxies with their counterparts in the low-density 
field, and its relation to the present study. Finally, we dis- 
cuss recent predictions for cluster-centric age gradients from 
N-body simulations and semi-analytic modelling. 

5.1 Cluster-centric gradients in stellar 
populations 

Our present study is, of course, not the first to address the 
issue of environmental effects on early-type galaxy proper- 
ties. In comparing to previous results, we concentrate on 
primarily spectroscopic tests, which can more easily sepa- 
rate age and metallicity effects. Among published studies, 
the closest in spirit to our approach are those of Guzman et 
al. (1992) and Carter et al. (2002). Each of these works used 
line-strengths for galaxies in only a single cluster, Coma, in 
contrast to our results which implicitly average over many 
clusters. 

Guzman et al. (1992) reported a zero-point offset of 

0. 017 ±0.005 mag in the Mg2 — a relation between a sample 
from the core of Coma and a 'halo' sample at larger radius. 
Their study included only morphologically-classified ellipti- 
cals, probing typically typically higher mass systems than 
our sample (90% with a > 150kms _1 , the median value for 
NFPS). The sense of the offset is the same as observed here, 

1. e. that at fixed velocity dispersion, the Mg line-strength is 
weaker at larger radius. With only one index, Guzman et 
al. could not distinguish among alternative interpretations 
of the offset as attempted here, but instead described the 
effect in terms only of age, with more recent star-formation 
occurring in the outer galaxies. The median radius for the 
two samples corresponds to 0.3i?2oo and 1.5i?200. Over this 
range, the radial term in our fit for the Mgb5177 index 
would imply a change of 0.18 A, equivalent to 0.011 mag 
in Mg2. Quantitatively, then, the radial behaviour of Mg 
found in this paper is marginally weaker than that claimed 
by Guzman et al. 

Carter et al. (2002) investigated residual gradients us- 
ing Mg2, Hbeta and (Fe) indices for a sample extending to 
1 R200 (80 arcmin) in Coma, without morphological or colour 
selection. They report a very significant radial trend in Mg2, 
with a change of 0.03 mag from core to 7?2oo- Carter et al. in- 
terpret their results as metallicity differences at fixed mass, 
and comment on the possibility of an abundance ratio ef- 
fect, since their measured Fe trend is weak. Converting to 
Mgb5177, the implied radial trend is ~ 0.6 ±0.2 A per -R200. 
Thus the Carter et al. trend, although in the same direc- 
tion as our result, is four times steeper (and three times 
steeper than Guzman et al. 1992) . Their Hbeta radial trend 
is also steeper than reported here. An important caveat to 
these comparisons is that the luminosity range spanned by 
the Carter et al. sample is larger than that used here, or by 
Guzman et al. We have already commented that the radial 
effects could be larger for galaxies of lower velocity disper- 
sion, which would work in the right direction to reconcile 
the two studies. 

Pimbblet et al. (2006) recently measured a radial shift 
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Figure 13. Predicted versus observed radial slopes for the restricted models, i.e. where only one or two of the stellar population 
parameters are allowed to vary with radius. Each panel is equivalent to Figure but for the best-fitting models of the class indicated. 
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in the B — R colour of the red sequence, in 11 clusters at 
z ~ 0.1. The fitted sample spans a luminosity range com- 
parable to NFPS. The magnitude of the colour trend is 
—0.011 ± 0.003 mag per Mpc. Given also their measure- 
ment of a significant rising radial trend in the H<5 absorption, 
Pimbblet et al. argue this is due to younger populations at 
larger radius. Since their definition of RS is not the same as 
the HdF index employed here, we compare on the basis of the 
B—R colour gradient predicted from our age, metallicity and 
[a/Fe] trends. Using models by Maraston (1998, 2005), our 
radial age gradient would yield S(B - R) = -0.029 ± 0.008 
blueing at the virial radius, while the metallicity effect con- 
tributes S(B-R) = +0.003±0.006 reddening. The Maraston 
models do not account for [a/Fe] variations; using a differ- 
ential estimate of g^f /^ ~ 0.12 derived from Salasnich et 
al. (2000), the reduced [a/Fe] at large radius further reddens 
the outer galaxies by +0.006 ± 0.001. The net effect is thus 
8(B -R) = -0.020 ± 0.010 at the virial radius. The median 
-R200 for the Pimbblet et al. sample is 2.9 Mpc. Thus our 
prediction for their colour trend is —0.007 ± 0.003 mag per 
h§Q Mpc, which is consistent with their result. 

In summary, our results from NFPS strongly confirm 
previous hints of radial differences in stellar populations, 
while our more general approach helps to resolve apparently 
contradictory claims concerning the driving parameters of 
the trends. 



5.2 Comparison with field vs cluster studies 

Arguably a more natural way to search for the influence 
of high-density environments is to compare galaxies in rich 
clusters to those in the 'field'. The caveat to this, of course, 
is that many alternative definitions of 'field' are in use, e.g. 
large distance to nearest neighbour; found in a large-scale 
under-density; not in a known rich cluster. As an example, 
an elliptical member of a compact group could qualify under 
some of these criteria, but not under others. 

It is tempting to compare the radial trends of galaxy 
populations in clusters with the results from field surveys, 
since galaxies falling into the cluster for the first time will 
not yet have experienced any effect of the cluster, and should 
match the field population. However, this comparison should 
be treated with some caution, since many galaxies may have 
been accreted into their host cluster as members of infalling 
groups, rather than as true field objects. If so, the outer 
cluster members may already have been substantially 'pre- 
processed' in the group environment (e.g. Kodama et al. 
2001), and should not be expected to match the field sam- 
ples. 

Many works have employed spectroscopic indices to 
test for systematic differences in the stellar populations of 
early-type galaxies in the field and in clusters. For example, 
Bernardi et al. (1998) found a small offset towards lower Mg2 
in the field, by 0.007±0.002 mag at given velocity dispersion, 
and interpreted this as an age difference with the field galax- 
ies 1 Gyr younger than cluster populations. Kuntschner et 
al. (2002) used multiple indices, reporting younger ages, by 
2-3 Gyr compared to clusters, for a small sample of very iso- 
lated galaxies. Thomas et al. (2005) obtained a similar offset 
of ~2 Gyr between field and cluster ellipticals. None of the 
above studies found evidence for an environmental shift in 



[a/Fe], although both Kuntschner et al. and Thomas et al. 
report an increased total metallicity [Z/H], which is not re- 
produced in our results. 

Most recently, Bernardi et al. (2006) analysed com- 
posite spectra derived from the Sloan Digital Sky Survey, 
for early-type galaxies selected by morphology and spectral 
type. They conclude that early-type galaxies in the lowest- 
density environments are younger by ~1 Gyr, than their 
counterparts of similar mass in the densest regions. The 
study finds a shift of 0.02 dex in [a/Fe], with relatively less 
a-enhancement in the field, but no significant offset is found 
for the overall metallicity. These results are qualitatively in 
striking agreement with the conclusions of this paper, if we 
allow that our outer cluster galaxies should approximately 
resemble their low-density population. 

5.3 Comparison with galaxy formation models 

The observed pattern of cluster-centric gradients agrees 
qualitatively with a model in which galaxies in the outer 
parts of cluster had more extended star formation activity, 
leading both to younger ages and more pollution of the inter- 
stellar gas by iron from SNe la, before star formation ceased. 
Thomas et al. (2005) constructed simple models of chemical 
enrichment in an extended burst, where the star-formation 
rate is Gaussian with FWHM r in Gyr, to describe the re- 
lationship between [a/Fe] and the formation time-scale r. 
Adopting their relation, [a/Fe]^ -| — glogr, our shift of 
A[a/Fe]= 0.05 between the core and the virial radius cor- 
responds to Alogr = 6 x A[a/Fe] « 0.3. Thus, galaxies at 
^200 formed their stars on time-scales twice as long as those 
for the inner galaxies. 

In high-resolution dark-matter simulations of cluster 
formation, sub-haloes found near the cluster centre at the 
present day were accreted at larger redshift. For instance 
Gao et al. (2004) find average infall redshifts of z ~ 0.4 at the 
virial radius, compared to z « 1 in the core. Together with 
the expectation that star-formation in the sub-haloes may 
be truncated after accretion to the cluster, perhaps through 
removal or exhaustion of the gas reservoir (Larson et al. 
1980; Balogh et al. 2000), this correlation will lead to an 
age gradient in the observed sense. Given a longer period of 
continuous star-formation prior to accretion, outer galaxies 
will also tend to have less a-enriched stellar populations. 

As a toy model, we consider at face value the average in- 
fall redshifts of Gao et al.'s Figure 15 (lower- left panel), and 
convolve with a star-formation rate which is constant until 
accretion, then declines exponentially with an e-folding time 
of 1 Gyr. The model yields present day luminosity-weighted 
ages of 10 Gyr in the centre, and 7 Gyr at the (physical) virial 
radius. Clearly a number of important factors have been ne- 
glected here: First, the average infall redshift shown by Gao 
ct al. gives the time since the sub-haloes entered the ulti- 
mate parent halo. In practice, some of these will have been 
absorbed within smaller parent haloes at earlier epochs, so 
that truncation may have taken place at a higher redshift. 
Second, as emphasized by Kodama & Bower (2001), the evo- 
lution of the global star-formation rate implies that the typ- 
ical field accreted galaxy at z « 0.4 was less active than that 
accreted atzssl (Madau, Pozzetti & Dickinson 1998). Fi- 
nally, and more trivially, the predicted gradients are in terms 
of physical radius from the cluster centre, rather than pro- 
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jected radius as observed. Each of these three effects leads 
to an over-estimation of the trends, so that the 30% age 
change from core to virial radius should be regarded as an 
upper limit on the likely gradient. 

Cosmological N-body simulations and semi-analytic 
galaxy formation models can provide more realistic descrip- 
tions of both the star-formation history and the merger his- 
tory of the accreted galaxies. With the large mass resolution 
afforded by the Virgo Consortium's Millennium Simulation 
(Springel et al. 2005), it is possible to follow sub-haloes cor- 
responding to sub-L* cluster members, and analyse their 
properties with reference to their location in the parent clus- 
ter. Recently, De Lucia et al. (2006) applied updated 'Mu- 
nich' semi-analytic recipes to Millennium Simulation merger 
trees, to model the formation of elliptical galaxies. They pre- 
dict that the luminosity-weighted stellar age declines from 
~12 Gyr at the cluster centres to ~10.5 Gyr at -R200, a ~13% 
effect (their Figure 8a). In comparing to our observational 
results, two caveats should be borne in mind. First, we mea- 
sure a trend in projected radius, whereas De Lucia et al. 
plot the physical distance. Second, De Lucia et al. do not 
correct for the overall age-versus-mass trend, and comment 
that some of their cluster-centric gradient arises from mass 
segregation. Examining their Figures 8a and 8c, most of the 
change in age is from 0.1 -R200 to 1.0 -R200, a range over which 
the median stellar mass of ellipticals changes by only ~2%. 
Even considering the inner O.I-R200, the 10% mass change, 
combined with their weak age-mass trend for massive galax- 
ies (their Figure 6a), cannot be primarily responsible for 
their radial gradient. Thus, despite the difficulties in match- 
ing the observed results to the models, our measured age 
shift of 15 ± 4% appears to be compatible with these predic- 
tions. In future work, we intend to make a more detailed con- 
frontation with the newly-updated 'Durham' semi-analytic 
predictions as described by Bower et al. (2005). 

In summary, a modest age-radius (and [a/Fe]-radius) 
relation is expected in hierarchical models, given the am- 
ple observational evidence for a quenching of star formation 
in the infall regions of nearby clusters (Lewis et al. 2002; 
Gomez et al. 2003). Our detection of a radial age trend in 
the red sequence, of the expected sense and order of mag- 
nitude, suggests that present-day red-sequence galaxies are 
the descendents of similar quenching mechanisms occurring 
at earlier cosmic epochs. 



6 CONCLUSIONS 

We have presented an analysis of the absorption line- 
strength indices of red-sequence galaxies in clusters, fo- 
cussing on the subtle effects of environment, as indicated 
by distance from the host cluster centre. 

Our principal conclusions can be summarised as follows: 

(i) The line-strengths of red-sequence galaxies exhibit 
systematic trends as a function of distance from the cen- 
tre of their host cluster, after controlling for the dominant 
trend with velocity dispersion. 

(ii) The cluster-centric trends are small compared to the 
overall changes of line-strengths with velocity dispersion, 
and small compared to the total scatter around these re- 
lations. None the less, for some indices the radial effects are 
highly significant. 



(iii) The pattern of gradients is different from that of the 
index versus velocity dispersion relations. The radial trends 
are not a scaled version of the mass trends. Whatever the 
underlying process driving the radial differences, it is dis- 
tinct from that which imposes the overall mass sequence. 

(iv) There is a segregation in the indices between those 
which respond strongly to age or [a/Fe] ratio, where the 
radial trends are opposite in sign to the mass trends, and 
the [Z/H]-indicators, which have radial trends of the same 
sign as the mass trends. 

(v) Within the context of the TMBK stellar population 
models, the cluster-centric gradients are best described by 
a significant variation in age and relative a-abundance. The 
sense is that red-sequence galaxies at one virial radius, -R200, 
on average have younger stars by 15 ±4% than red-sequence 
galaxies of the same mass at the cluster centre. The a abun- 
dance, relative to iron, is 10 ± 2% lower at -R200 than at the 
cluster centre, again at fixed mass. The overall metallicity 
[Z/H] does not show a significant trend. 

(vi) Available morphology data suggests that although 
stellar populations appear to depend strongly on bulge- 
fraction B/T, this does not translate into a significant con- 
tribution to the radial trend. The radial behaviour is not 
driven by a changing morphological mix, but rather a change 
in typical stellar populations for galaxies of the same mor- 
phology, as well as the same mass. 

(vii) Robustness tests using subsets of the input data 
suggest that the radial trend signals are quite evenly di- 
vided between distant and nearby clusters, and between the 
two celestial hemispheres. Using only the relaxed clusters 
with well-defined centres, we recover slightly stronger trends. 
There is a suggestion that the radial trend in age is strongest 
for galaxies of lower mass. 

(viii) The age trend determined here from the observed 
line-strength indices is in the sense generically predicted 
from hierarchical clustering with star-formation truncation 
in clusters. Numerically, there is good agreement with very 
recent predictions from semi-analytic galaxy formation mod- 
els (De Lucia et al. 2006), which predict an age decline of 

~13% OUt to i?200- 



In summary, we have shown that the average star for- 
mation histories of red-sequence galaxies depend on their lo- 
cation in the cluster, in addition to the overall dependence 
on galaxy mass. Our results lend support to the notion that 
the red sequence was established through the accretion of 
star-forming galaxies from the field, with a subsequent envi- 
ronmental truncation of star formation. Together with our 
earlier measurement of the mass dependence of the aver- 
age star- formation history (Nelan et al. 2005), a picture is 
beginning to emerge for how the red sequence became popu- 
lated as a cosmic history progressed. Further progress in this 
area will be driven by improved spectroscopic observations 
of faint cluster members, which appear to exhibit stronger 
signatures of later accretion. Additional insights should re- 
sult from comparing in greater detail with morphological 
and photometric evidence for galaxy transformation in clus- 
ters (e.g. Stevenson et al. 2006), with observations at higher 
redshifts (e.g. Tanaka et al. 2005), and with improved mod- 
elling of the hierarchical galaxy formation process. 



20 Russell J. Smith et al. 



ACKNOWLEDGMENTS 

RJS was supported for part of this work from the PPARC 
rolling grant 'Extragalactic Astronomy & Cosmology at 
Durham 2005-2010' (C.S. Frenk). MJH acknowledges sup- 
port from the NSERC of Canada. JEN was supported by 
a Dartmouth College Fellowship, a Graduate Assistance in 
Areas of National Need fellowship, and a NASA space grant. 
We gratefully acknowledge the substantial assignment of 
NOAO observing resources to the NFPS program. This pub- 
lication makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analy- 
sis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the Na- 
tional Science Foundation. 



REFERENCES 

Balogh M.L., Navarro J.F., Morris S.L. 2000, ApJ, 540, 113 
Balogh M.L., Schade D., Morris S.L., Yee H.K.C., Carlberg R.G., 

Ellingson E. 1998, ApJ, 504, L75 
Bernardi M., Renzini A., da Costa L.N., Wegner G., Alonso M.V., 

Pellegrini P.S., Rite C, Willmer C.N.A. 1998, ApJ, 508, L143 
Bernardi M., Nichol R.C., Sheth R.K., Miller C.J., Brinkmann J. 

2006, AJ, 131, 1288 
Bower R.G., Benson A.J., Malbon R., Helly J.C., Fren k C.S., 

Baugh CM., Cole S., Lacey C.G., astro-ph/0511338 
Caldwell N., Rose J.A., Concannon K.D. 2003, AJ, 125, 2891 
Carlberg R.G. et al. 1997, ApJ, 485, L13 
Carter D. et al. 2002, ApJ, 567, 772 
Cole S.M. et al. 2001, MNRAS, 326, 255 

Crawford C. S., Allen S. W., Ebeling H., Edge A. C, Fabian A. 

C. 1999, MNRAS, 857 
De Lucia G., Springel V., White S.D.M., Croton D., Kauffman 

G. 2006, MNRAS, 366, 499 
Djorgovski S., Davis M. 1987, ApJ, 313, 59 
Dressier A. 1980, ApJ, 236, 772 

Gao L., White S.D.M., Jenkins A., Stoehr F., Springel V. 2004, 

MNRAS, 355, 819 
Guzman R., Lucey J.R., Carter D., Terlevich R.J. 1992, MNRAS, 

257, 187 

Gomez, P. et al. 2003, ApJ, 584, 210 

Gonzalez J.J. 1993, Ph.D. thesis, Univ. California, Santa Cruz 
Goudfrooij P. 1999, ASP Conf. Ser. 163: Star Formation in Early 

Type Galaxies, 163, 55 
Hogg DW. et al. 2002, AJ, 124, 646 
Hogg DW. et al. 2004, ApJ, 601, L29 
Hubble E., Humason, M.L. 1931, ApJ, 71, 43 
Kodama T., Bower R.G. 2001, MNRAS, 321, 18 
Kodama T., Smail I., Nakata F., Okamura S., Bower R.G. 2001, 

ApJ, 562, L9 

Jarrett T. H., Chester T., Cutri R., Schneider S., Skrutskie M., 

Huchra J. P. 2000, AJ, 119, 2498 
J0rgensen I., Franx M., Kjasrgaard P. 1995, MNRAS, 276, 1341 
Kochanek C.S., Pahre M.A., Falco E.E., Huchra J. P., Mader J., 

Jarret T.H., Chester T., Cutri R., Schneider S.E. 2001, ApJ, 

560, 566 

Kuntschner H., Lucey J.R., Smith R.J., Hudson M.J., Davics R.L. 

2001, MNRAS, 323, 615 
Kuntschner H., Smith R.J., Colless M.M., Davies R.L., Kaldare 

R., Vazdekis A. 2002 MNRAS, 337, 172 
Larson R.B., Tinsley B.M., Caldwell C.N. 1980, ApJ, 237, 692 
Lee H.-c, Worthey G. 2005, ApL160, 176 
Lewis, I. et al. 2002, MNRAS, 334, 673 
Madau P., Pozzetti L., Dickinson M. 1998, ApJ, 498, 106 



Maraston C. 1998, MNRAS, 300, 872 
Maraston C. 2005, MNRAS, 362, 799 

Mehlert D., Thomas D., Saglia R. P., Bender R., Wegner G., 2003 

A&A, 407, 423 
Melnick J., Sargent W.L.W. 1977, ApJ, 215, 401 
Moore B., Diemand J., Stadel J. 2004, IAU Colloq. 195: Outskirts 

of Galaxy Clusters: Intense Life in the Suburbs, 513 
Nelan J.E., Smith R.J., Hudson M.J., Wegner G.A., Lucey J.R., 

Moore SAW., Quinney S.J., Suntzeff N.B. 2005, ApJ, 632, 

137 

Pimbblet K.A., Smail I., Edge A.C., O'Hely E., Couch W.J., 

Zabludoff A.I. 2006, MNRAS, 366, 645 
Poggianti B.M., et al. 2001, ApJ, 562, 689 
Proctor R.N., Samson A.E. 2002, MNRAS, 333, 517 
Proctor R., Forbes DA., Hau G.K.T., Beasley M.A., De Silva 

G.M., Contreras R., Terlevich A.I. 2004, MNRAS, 349, 1381 
Quint ero A.D., Berli nd A., Blanton M.R., Hogg DW. 2005, 

astro-ph/0512004 
Salasnich B., Girardi L., \Veiss A.. Chiosi C. 2000. A&A, 361, 

1023 

Smith R.J. 2005, MNRAS, 359, 975 

Smith R.J. et al. 2004, AJ, 128, 1558 

Springel V. et al., 2005, Nature, 4435, 629 

Stevenson J.B., Hudson M.J., et al. 2006, in preparation 

Tanaka M., Kodama T., Arimoto N, Okamura S., Umetsu K., 

Shimasaku K., Tanaka I., Yamada T. 2005, MNRAS, 362, 

268 

Terlevich A.I., Caldwell N. & Bower R.G. 2001, MNRAS, 326, 
1547 

Thomas D., Davies R.L. 2006, MNRAS, 366, 510 
Thomas D., Maraston C, Bender, R. 2002, MNRAS, 339, 897 
Thomas D., Maraston C, Korn, A. 2004, MNRAS, 351, L19 
Thomas D., Maraston C, Bender R., Mendes de Oliveira C. 2005, 
ApJ, 621, 673 

Trager S.C., Worthey G., Faber S.M., Burstein D., Gonzalez J.J. 

1998, ApL, 116, 1 
Trager S.C., Faber S.M., Worthey G., Gonzalez J.J. 2000a, A J, 

119, 1645 

Tripicco M.J., Bell RA. 1995, AJ, 110, 3035 
Visvanathan N., Sandage A. 1977, ApJ, 216, 214 
Worthey G. 1994, ApL95, 107 
Worthey C, Ottoviani D.L. 1997, ApLlll, 377 
Worthey G., Faber S. M., Gonzlez J. J., 1992, ApJ, 398, 69 
Yan R., Newman J. A., Faber S.M., Konidaris N., Koo D., Davis 
M., 2005 |astro-ph/0512446| 



